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Abstract: 

We investigate structural, kinematical, and chemical properties of stars and gas in the Small Magellanic 
Cloud (SMC) interacting with the Large Magellanic Cloud (LMC) and the Galaxy based on a series 
of self-consistent chemodynamical simulations. We adopt a new "dwarf spheroidal model" in which 
the SMC initially has both old stars with a spherical spatial distribution and an extended HI gas 
disk. We mainly investigate SMC's evolution for the last ~ 3 Gyr within which the Magellanic stream 
(MS) and the Magellanic bridge (MB) can be formed as a result of the LMC-SMC-Galaxy interaction. 
Our principal results, which can be tested against observations, are as follows. The final spatial 
distribution of the old stars projected onto the sky is spherical even after the strong LMC-SMC-Galaxy 
interaction, whereas that of the new ones is significantly flattened and appears to form a bar structure. 
Old stars have the line-of-sight velocity dispersion (ct) of ~ 30 km s~^ and slow rotation with the 
maximum rotational velocity (V^) of less than 20 km s^^ and show asymmetry in the radial profiles. 
New stars have a smaller a than old ones and a significant amount of rotation {V/a > 1). HI gas 
shows velocity dispersions of a = 10 — 40 km s~^, a high maximum rotational velocity {V ^ 50 km 
s~^), and the spatial distribution similar to that of new stars. The new stars with ages younger than 
3 Gyr show a negative metallicity gradient in the sense that more metal-rich stars are located in the 
inner regions of the SMC. The MB inevitably contains old stars with the surface mass densities of 
6 — 300 X IO^Mq deg~^ depending on initial stellar distributions of the modeled SMC. We find that 
the dwarf spheroidal model can explain more self-consistently the observed kinematical properties 
of stars and gas, compared with another type of the model ("the disk model") in which the SMC 
initially consists of stellar and gas disks. We suggest that the SMC needs to be modeled as hav- 
ing a spheroidal component rather than being a purely disk for better understanding the SMC evolution. 
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1 Introduction 

Dynamical interaction between the SMC, the LMC, 
and the Galaxy have long been considered to be closely 
associated not only with the evolution of the SMC 
(e.g., Yoshizawa & Noguchi 2003; YN03) but also with 
the formation of the MS (e.g., Putman et al. 1998) and 
the MB (e.g., MuUer et al. 2004). The observed struc- 
tural, kinematical, and chemical properties of stars, 
star clusters, and gas in the SMC have been exten- 
sively discussed in variously different aspects of the 
SMC, such as the star formation history (e.g., Rus- 
sell & Dopita 1992; Pagel & Tautvaisiene 1998; Pi- 
atti et al. 2007), the total dynamical mass (e.g., Do- 
pita et al. 1985; Stanimirovic et al. 2004, S04), and 
the interaction history with the LMC (YN03; Bekki & 
Chiba 2007). One of the key, longstanding problems 
related to the SMC as a galaxy is to understand how 
the past dynamical interaction between the SMC, the 
LMC, and the galaxy could have influenced the chem- 
ical and dynamical evolution and the star formation 
history of the SMC (e.g., Westerlund 1997). 

In order to solve the above mentioned key problem. 



we need to know not only three-dimensional (3D) space 
motions of the Magellanic Clouds (MCs) but also dy- 
namical and chemical properties of stars with different 
ages and metallicities. Recent proper motion measure- 
ments of the MCs and numerical studies of the orbits 
of the MCs based on the proper motion measurements 
have suggested that the MCs can be a unbound with 
each other and they are unbound to the Galaxy (Kalli- 
vayaUl et al. 2006, K06; Besla et al. 2007). The two- 
dimensional (2D) distributions of stellar populations 
with different ages and metallicities derived from AGB 
stars have revealed star formation histories in different 
local regions of the SMC (Cioni et al. 2006). Observa- 
tional studies of stellar populations of the SMC based 
on their color-magnitude diagrams have revealed not 
only the 2D structures dependent on ages and metal- 
licities but also a more precise global star formation 
history of the SMC (e.g., Harris & Zaritsky 2004). 
Thus, these observations on 2D structures of different 
stellar populations, if compared with the correspond- 
ing theoretical studies, can provide a vital clue to the 
above-mentioned key question of the SMC. 

Previous theoretical models, however, did not pro- 
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vide useful predictions on the 2D (or 3D) structures 
of stars with different ages and metallicities owing to 
the strong limitations in their numerical simulations. 
YN03 discussed the star formation history of the SMC 
and the origin of structural and kinematicaJ proper- 
ties of the MS. They however did not discuss chemical 
and dynamical properties of stars and gas within the 
SMC. Although previous one-zone chemical evolution 
models discussed the long-term chemical evolution his- 
tory of the SMC (e.g., Spite et al. 1988; Tsujimoto et 
al. 1995; Idiart et al. 2007), they did not discuss 
dynamical properties of the SMC. Thus no theoret- 
ical studies have been done which attempt to inves- 
tigate both dynamical and chemical properties of the 
SMC and thereby provide reasonable explanations for 
the above-mentioned observations. Given that only 
chemodynamical simulations can provide predictions 
on 3D-dynamical and chemical properties in a self- 
consistent manner, it is vital for theoretical studies 
to investigate the above-mentioned key observational 
properties based on chemodynamical simulations of 
the SMC. 

The purpose of this paper is thus to numerically 
investigate physical properties of the SMC in a fully 
self-consistent manner and thereby try to provide an- 
swers for the above-mentioned unresolved problems. 
We mainly discuss the origin of structural, kinemati- 
cal, and chemical properties of the SMC's stellar and 
gaseous components derived in previous observations 
(e.g., Suntzeff et al. 1986; Torres & Caranza 1987; 
Hardy et al. 1989; Russell & Dopita 1992; Hatzidim- 
itiriou et al. 1993; S04). We also discuss the latest 
observations that previous theoretical works did not 
discuss: the structure of the stellar halo (Noel & Gal- 
lart 2007) , kinematical properties of older stellar pop- 
ulations (e.g., Harris & Zarisky 2006), possibly purely 
gaseous streams in the MB (Harris 2007), and bifur- 
cated structures and kinematics of the MC (Matthews 
& Staveley-Smith 2007). We however do not intend to 
discuss some important aspects of the SMC evolution, 
such as the long-term evolution of the MCs and the 
formation of globular cluster (GCs), which were dis- 
cussed to some extent in our previous papers (Bekki 
et al. 2004; Bekki & Chiba 2005; paper I, Bekki 2007) 
in the preset study. We plan to discuss these in our 
paper HI (Bekki & Chiba 2008a) based on more so- 
phisticated and self-consistent models not only for the 
MCs but also for the Galaxy. 

The layout of this paper is as follows. In §2, we 
summarize our numerical models used in the present 
study and describe the methods for analyzing struc- 
ture and kinematics of the simulated SMC. In §3, we 
present numerical results on the time evolution of mor- 
phology, metallicity distribution, and dynamical prop- 
erties of the SMC. In §4, we discuss the above-mentioned 
outstanding issues related to formation and evolution 
of the SMC. The conclusions of the preset study are 
given in §5. Since the origin of the MS has been al- 
ready discussed extensively in previous numerical sim- 
ulations (e.g., Gardiner & Noguchi 1996, GN; Con- 
nors et al. 2006, C06), we briefly discuss it in the 
present paper. The present paper is still based on sim- 
ilar approaches used in previous numerical works (e.g.. 



simple analytic models for dynamical friction and the 
fixed isothermal density profile for the Galactic dark 
matter halo) . The more sophisticated models in which 
the LMC, the SMC, and the Galaxy are represented 
by fully consistent N-body models will be described in 
our forthcoming papers (Bekki & Chiba 2008a). 
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Figure 1: Orbital evolution of tlie MCs for the 
models CI and C2 with the GN (upper) and the 
K06 (lower) velocity types. The time evolution 
of distances between MCs (thick solid), the LMC 
and the Galaxy (thin solid), and the SMC and 
the Galaxy (dotted) are separately shown. Here 
Mlmc = 2.0 X lOi°M0 and Msmc = 3.0 x IO^Mq 
are assumed for the two models. 



2 Model 

The present investigation is two-fold. First, we in- 
vestigate orbital evolution of the MCs with respect to 
the Galaxy for a given set of initial velocities of the 
MCs by using a backward integration scheme (Murai 
& P\ijimoto 1980; MF). Then we investigate chemo- 
dynamical evolution of the SMC on the orbits by us- 
ing our original chemodynamical simulations with star 
formation models. Since the details of the methods 
for orbital calculations and numerical techniques for 
chemodynamical simulations of the MCs are given in 
the paper I, wc briefly describe these in the present 
paper. The present model is significantly better than 
the previous ones (e.g., GN and C06) in the sense that 
(i) the SMC model is much more consistent with the 
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observations and (ii) star formation and chemical evo- 
lution processes are better modeled. 

2.1 Derivation of the orbits of the 
MCs 

The gravitational potential of the Galaxy $g is as- 
sumed to have the logarithmic potential; 

3.G(r) = -yo'lnr, (1) 

whore r and Vo arc the distance from the Galactic cen- 
ter and the constant rotational velocity (= 220 km 
s~^), respectively. We consider this potential, firstly 
because the present model can be consistent with that 
of the paper I (and those of previous papers, e.g., MF) 
and secondly because we can estimate the drag force 
due to dynamical friction between the MCs and the 
Galaxy in an analytical way for the orbital calculations 
of the MCs. Thanks to this adoption of the potential, 
we can compare the present results with previous ones 
(e.g., GN, C06, K06, and Besla et al. 2007). 

Owing to the adopted isothermal density profile of 
the halo (i.e., ever increasing mass with radius), the 
MCs are more likely to show bound orbits for given 
initial velocities. If we adopt that the latest model 
for the Galaxy with a more realistic halo model (e.g., 
Widrow & Dubinski 2005), the orbital evolution of the 
MCs can be different owing to the smaller mass of the 
Galaxy in its outer part {R > 100 kpc). However, it is 
not practical for the present study to adopt the model 
by Widrow & Dubinski (2005), because analytical es- 
timation of the drag force due to dynamical friction 
between the MCs and the Galaxy is not feasible. 

The LMC is assumed to have the Plummer poten- 
tial; 

$L(rL) = -MLMc/(rL' + alf\ (2) 

where Mlmc, '"l, and aL arc the total mass of the 
LMC, the distance from the LMC, and the cH'ective 
radius, respectively. We adopt the same value of ttL 
(=3 kpc) as previous numerical studies adopted (e.g., 
GN). The SMC is also assumed to have the Plummer 
potential; 

Mrs) = -GMsMc/(rs' + alf\, (3) 

where Msmc , ^"3 , and as are the total mass of the SMC, 
the distance from the SMC, and the eflFective radius, 
respectively. We adopt the same value of as {—2 kpc) 
as previous numerical studies adopted (e.g., GN). 

We consider the dynamical friction due to the pres- 
ence of the Galactic dark matter halo both for the 
LMC-Galaxy interaction and for the SMC-Galaxy one 
and adopt the following expression (Binney & Tremaine 
1987); 

FfricG = -0.428 In Ag , (4) 

where r is the distance of the LMC (the SMC) from 

the center of the Galaxy. The mass AI is either A^lmc 
or AfsMC, depending on which Cloud's orbit (i.e., LMC 
or SMC) we calculate. We adopt the reasonable value 
of 3.0 for the Coulomb logarithm Ag (GSF; GN) both 



in the orbital calculation of the LMC and in that of 
the SMC. 

By integrating equations of the motions of the MCs 
toward the past from the present epoch, we investigate 
orbital evolution of the MCs for given Mlmc, AfsMC, 
initial positions and velocities of the MC. We adopt 
the reasonable sets of orbital parameters that are con- 
sistent with observations. The current Galactic coor- 
dinate (6, Z), where I and b are the Galactic longitude 
and latitude, respectively, is (—32.89, 280.46) for the 
LMC and (-44.30, 302.79) for the SMC. Therefore the 
current positions {X,Y,Z) of the LMC, {Xl,Yl,Zl) 
in units of kpc, are (—1.0, —40.8, —26.8) and those of 
the SMC, (Xs,Ys,Zs), are (13.6,-34.3,-39.8). The 
current distance and the Galactocentric radial velocity 
of the LMC (SMC) is 80 (7) km s'^ 

The current space velocities or {U, V, W) in units 
of km s~^ are the most important parameters that de- 
termine the orbital evolution of the MCs in the present 
models. They are represented by {U-l,Vl,Wl) for the 
LMC and by {Us,Vs,Ws) for the SMC. We mainly 
investigate the following two velocity types: The "GN 
type" with (C7l, Fl, Wl) = (-5, -225, 194) and {Us, Vs, Ws) 
= (40, -185, 171) and the "K06 type" with {Ul, Vl,W^) 
= (-86, -268, 252) and {Us,Vs, Ws) = (-87, -247, 149). 
The former is the same as those used in the models by 
GN whereas the latter is consistent with observations 
by K06. 

Fig. 1 describes the past ~ 3 Gyr orbital evolu- 
tion of the MCs in the models with the GN and K06 

velocity types for which Mlmc ~ 2 x 1O"'^''M0 and 
Msmc = 3 x 10'"'Mq are assumed. Here negative val- 
ues of the time, T, represent the past, with T = 
corresponding to the present epoch. As shown in this 
figure, the present orbital period of the MCs about the 
Galaxy is ~ 1.5 Gyr for the adopted gravitational po- 
tential and the masses of the MCs. The strong LMC- 
SMC interaction around 1.5 Gyr can be seen only in 
the model with the GN type, which suggests that the 
MS can not be formed from the SMC in the model 
with the K06 type. 

For a given set of initial positions and velocities 
of the MCs predicted by the above backward integra- 
tion scheme, we dynamically evolve the MCs by using 
GRAPE-SPH simulations (i.e., integrating the equa- 
tion of motion forward) for a fixed potential of the 
Galaxy. Since the mass of the SMC is very small in 
comparison with the Galaxy, we consider that the tidal 
effects of the Galaxy on the SMC (e.g., stripping of 
stars from the SMC) can be properly treated by our 
models. Dynamical evolution of the SMC under the 
live potential of the Galaxy will be investigated by our 
future studies and the results will be compared with 
those of the present study. 

We mainly discuss the models with "classical" bound 
orbits of the LMC and the SMC (i.e., GN-type) around 
the Galaxy in order to reproduce well the observed 
properties of the MS. The bound orbits are also sug- 
gested by Shattow & Loeb (2008) which adopted Vo a 
factor of 14% larger than 220 km s~^ adopted in the 
present study. They also suggested that the bound 
orbits can be consistent with the latest proper motion 
measurements by K05 as long as the higher circular ve- 
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locity is adopted. Thus our main investigation of the 
GN-type bound orbits can be justified in the present 
study. 
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Figure 2: Circular velocities (14) dependent on 
the distance (R) from the center of the SMC (up- 
per) and total masses within R (lower) for the 
model G5 with Mt = 3.0 x IO^Mq (a) and Gl 
with Mt = 4.5 X 1O^M0 (b). K is calculated sepa- 
rately for the mass profiles of the dark matter halo 
(solid), the old stars (dotted), the gas (short-dash), 
and all components (long-dash). 



2.2 GRAPE-SPH chemodynamical sim- 
ulations 



For a given orbit of the SMC, we investigate chemody- 
namical evolutions of the SMC by adopting the fol- 
lowing two different galaxy models: (i) the "dwarf 
spheroidal" model in which the SMC has a stellar spheroidal 
and an extended HI gas disk and (ii) the "disk" one 
in which the SMC has stellar and gaseous disks. Al- 
though the SMC is classified as an irregular dwarf of 
DDO type Ir IV- V (van den Bergh 2000), the usage of 
"dwarf" spheroidal for the SMC would not be so ap- 
propriate, given that the SMC's stellar mass is as large 
as IO^Mq (S04): we use this term just for convenience. 

Although the above disk model was adopted in pre- 
vious models of the SMC (MK, GN, YN, and C06), 
we mainly investigate the dwarf spheroidal model in 
the present study. This is mainly because the dwaxf 
spheroidal model can better explain the observed struc- 
tural and kincrnatical properties of the SMC. The LMC 
is represented by a test particle with the mass of Mlmc 
and the Plummer potential described in §2.1 for most 
models. The results of simulations in which both the 
LMC and the SMC are represented by gaseous and 
stellar particles in a fully self-consistent manner will 
be described in (Bekki & Chiba 2008a) with main 
points of analysis significantly different from those of 
the present one. 
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Figure 3: The final distribution of gas projected 
onto the y-z plane for the model CI. The orbits 
of the LMC (triangle) and the SMC (square) with 
respect to the center of the Galaxy are shown by 
solid and dotted lines, respectively. For compari- 
son, the Galaxy with the disk size of 17.5 kpc is 
shown by a thick solid line. 



2.2.1 The dwarf spheroidcd model 

A significant fraction of dwarf elliptical and spheroidal 
galaxies arc observed to have exponential luminosity 
profiles (e.g., Icliikawa, Wakamatsu, & Okamura 1986 
■Andredakis & Sanders 1994): we do not discriminate 
dwarf ellipticals and spheroids just for convenience in 
the present paper. Therefore the projected radial den- 
sity profile (Esph) of the stellar spheroid with the mass 
of Afsph and the scale length of Ogph is assumed to be 
an exponential one: 



Ssph(-R) = Ssph,o X exp(-i?,/agph). 



(5) 



where Esph,o and R are the central surface density and 
the projected radius from the center of the spheroid 
(i.e., the center of the SMC), respectively. In order to 
derive the three dimensional (3D) density field {psph{'r), 
where r is the distance from the center of the spheroid) 
from Esph, we can use the following formula (Binney 
& Tremaine 1986): 



Psph(r-) 



dS,ph(i?) dR 



dR VJi^ - r2 ■ 



(6) 



We rmmerically estimate the Pspii('") profile for a 
given Esph(i?) in a model with Afsph and asph- The 
outer cut-off radius of the stellar spheroid (-Rsmc,s) 
is assumed to be a free parameter. The spheroid is 
assumed to be embedded in a massive dark matter 
halo with the total mass of Afdm and the truncation 
radius of -Rsmc- This -Rsmc represents the size of the 
SMC within which almost all particles are included. 
The radial density profile is assumed to have either the 
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Table 1: Model parameters for N-body simulations 



Model name^ 


MsMC 1 


Velocity types ^ 


Morphology ^ 


-RsMC 


1 


f 9 
Js — 


J, h 




Comments 


Gl 


0.45 


GN 


dSph 


7.5 


0.5 


1.0 


2.0 


250 


fiducial 


G2 


0.45 


GN 


dSph 


7.5 


0.25 


1.0 


2.0 


250 


more compact 


G3 


0.45 


GN 


disk 


7.5 


0.25 


1.0 


4.0 


250 


disk model 


G4 


0.45 


GN 


dSph 


7.5 


0.5 


1.0 


2.0 


250 


2/met = 0.002 


G5 


0.30 


GN 


dSph 


5.0 


0.5 


2.0 


2.0 


250 


lower mass 


G6 


0.30 


GN 


dSph 


7.5 


0.25 


1.0 


2.0 


250 




G7 


0.45 


K06 


dSph 


7.5 


0.5 


1.0 


2.0 


250 


orbit from K06 


11 


0.45 




dSph 


7.5 


0.5 


1.0 


2.0 


250 


isolated 


12 


0.45 




dSph 


7.5 


0.5 


1.0 


2.0 


10^ 




13 


0.45 




dSph 


7.5 


0.5 


1.0 


2.0 


10* 




CI 


0.30 


GN 


dSph 


7.5 


0.25 


2.0 


2.0 


250 


coUisionless 


C2 


0.30 


K06 


dSph 


7.5 


0.25 


2.0 


2.0 


250 




C3 


0.45 


GN 


dSph 


7.5 


0.5 


1.0 


2.0 


250 




C4 


0.30 


GN 


dSph 


5.0 


0.5 


2.0 


2.0 


250 





""G" represents gaseous models in which gas physics, star formation in gas, and chemical enrichment processes are all included 
whereas "C" represents coUisionless models in which such processes are switched off. "I" represents isolated models in which interaction 
between the MCs and the Galaxy is not included at all. 

'The total mass of the SMC in units of 10^° Mq. 

cQN" and "K06" describes the orbits of the MCs that are consistent with those by GN and K06, respectively. 

■^The initial morphology of old stars: "dSph" and "disk" represent the dwarf spheroidal and the disk models, respectively. We use 
the term of dwarf spheroid for convenience, though it would not be so appropriate for the SMC to be classified as a "dwarf" spheroid 
for its luminosity. 

'The size of the SMC which corresponds to the tidal radius in units of kpc. 
^The ratio of the stellar disk size (-RsMCs) to the SMC one (-RsMc)- 
^The mass ratio of gas to stars in the SMC. 

''The size ratio of the gaseous disk to the stellar spheroid or disk in the SMC. 

'The initial gaseous temperate in units of K. The low temperature of cold gas is adopted so that the outer disk can keep its initial 
thin configuration for most models. 
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universal "NFW" one with the central cusp (Navarro, 
Frenk & White 1996) or a flat "core" (Burkcrt 1995; 
Salucci & Burkert 2000; referred to as the SB profile 
from now on). The NFW profile is described as: 



Pdm(r) 



Po 



(r/r,)(l+r/rs)2 



(7) 



where r, po, and rs are the distance from the center of 
the cluster, the scale density, and the scale-length of 
the dark halo, respectively. The SB profile is described 
as: 



Pdm(r-) = 



Pdmfi 



(8) 



(r-|-adm)(r2-|-0dm^)' 
where pdmfi and adm are the central dark matter den- 
sity and the core (scale) radius, respectively. The halo 
is truncated at 7.5 kpc from the center of the SMC. 

S04 demonstrated that (i) the HI rotation curve 
is very slowly rising within the central few kpc of the 
SMC with the maximum circular velocity of ~ 60 km 
s"^ and (ii) the total mass of the SMC is ~ 2.4x lO^M© 
within the central ~ 3 kpc. These observed properties 
of the SMC are better explained by the SB profile (see 
more detailed discussions on this for Bekki & Stan- 
imirovic 2008): the NFW proflle with the reasonable 
parameters for the SMC (e.g., c=5 — 20, whore c is the 
ratio of the virial radius to r^) is not so consistent with 
the above observation (i) . We find that the SB profiles 
with adm as large as 5 — 7 kpc can be more consistent 
with the observations for Msmc = 3 — 4.5 x lO^M© , 
though the slowly rising profile can not be exactly the 
same as the observed. We therefore show the results 
of the models with the SB profiles. We assume that 
Afdm/Afsph is 2.33 and the halo is truncated at the 
tidal radius (n) of the SMC. The previous models for 
the SMC assumed that n ranges from 5 kpc (GN) to 
7.5 kpc (C06). JisMC is the same as n in the present 
model and regarded as a free parameter. 

The stellar spheroid is assumed to be supported 
purely by velocity dispersion and its dispersion is as- 
sumed to be isotropic. We therefore estimate the ve- 
locities of old stellar particles from the gravitational 
potential at the positions where they are located. In 
detail, we first calculate the one-dimensional isotropic 
dispersion according to the following formula: 



U{r) 



(9) 



where U{r) is the gravitational potential at the posi- 
tion r which is determined by the mass distributions 
of the dark matter halo and the stellar spheroid. Then 
we allocate a velocity to stellar particle so that the 
distribution of velocities of these particles can have a 
Gaussian form with a dispersion equal to cr^(r). We 
adopt this method as an approximation and thus need 
to let the halo dynamically relaxed to obtain a real 
dynamically equilibrium for a given distribution of the 
SMC: this approximation is quite useful, in particu- 
lar, for implementing the triaxial potential of the disk 
model later described. 

The SMC is assumed to have an extended HI gas 
disk with the total gas mass of AfsMc.g, the size of 
RsMCs, and the exponential profile for the projected 
gaseous density. The scale length of the disk is set 



to be 0.2_RsMC,g for all spheroidal models. The size 
ratio (rg) of -RsMC.g to -Rsmc.s is a key parameter 
which determines whether the MS is composed almost 
purely of gas (YN03). Furthermore HI diameters of 
gas-rich galaxies axe generally observed to be signifi- 
cantly larger than their optical disks (Broeils & van 
Woerden 1994). Guided by these previous works, we 
investigate models with 2 < r-g < 4 in the present 
study. The mass ratio (/g) of MsMC.g to Msmc.s is 
a key parameter which determines the masses of the 
MS and the LA (YN03). We find that /g should be 
at least ~ 2.0 for Msmc = 3 x lO^M© and ~ 1.0 for 
AfsMC = 4.5 X 1O®M0 to explain the total gas mass in 
the MB and the MS regions. We thus show the results 
of the models with 1 < /g < 2 in the present study, 
though we have run models with /g < 0.5 as done in 
GN and C06. 

Since the disk component in this spheroidal model 
is composed purely of gas, it is highly likely to be dy- 
namically supported by rotation: gas dynamics (e.g., 
cooling and shock) in the disk can quickly dissipate 
away kinematic energy resulting from random motion 
of gas. Such a cold gas disk, however, can be quite 
unstable against dynamical instability (e.g., ring or 
bar instability) , if the mass fraction of the disk is sig- 
nificant (e.g., Binney & Tremairie 1987). For such a 
massive gas disk, the disk may well have a significant 
amount of kinetic energy due to random motion of 
the gas. In addition to the rotational velocity made 
by the gravitational field of the SMC, the initial ra- 
dial and azimuthal velocity dispersion are thus given 
to the gaseous component according to the epicyclic 
theory with Toomre's parameter (Binney & Tremaine 
1987) Q = 1.5: thus the cold gas is stabilized by its 
random motion. 

Since the SMC has a massive gaseous disk sur- 
rounding old stars, we need to adopt the following 
unique procedures to construct a dynamical equilib- 
rium of the SMC. Firstly we give velocities of dark 
matter particles for given density profiles of the dark 
matter halo and the stellar spheroid. Secondly we let 
the dark matter halo evolve for ten dynamical time 
scales with the spheroid being fixed. During this re- 
laxation process, the mass distribution of the halo can 
slightly change so that the system can reach a new 
dynamical equilibrium. Thirdly we allocate velocities 
for stellar particles according to the equation (9) and 
add a gas disk to the system. Fourthly, we again let 
the system evolve for ten dynamical time scales with 
the gas disk being fixed. During this relaxation pro- 
cess, the mass distribution of the spheroid (and the 
halo) can slightly change so that the system can reach 
a new dynamical equilibrium. Finally, we give rota- 
tional velocities to gaseous particles according to the 
circular velocities at their position for the new mass 
profile of the system. 

Fig. 2 shows the rotation curves of two dwarf 
spheroidal models with Msmc = 3.0 x lO^M© and 
i?SMC = 5.0 kpc and with Msmc = 4.5 x lO^M© and 
-RsMC ~ 7.5 kpc. These models show the maximum 
circular velocities of ^-^ 60 km s^^ , which is consistent 
with HI kinematics of the SMC (e.g., S04). For mod- 
els with small SMC masses (MsMC ~ 3.0 X lO'^Mo), 
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the initial SMC size should be as small as 5 kpc so 
that the final total mass of the SMC after efficient 
tidal stripping of stars and gas can be more consistent 
with the observations by S04. It should be stressed 
here that the rotation curves in this smaJl-mass model 
reach their peaks at ~ 1 kpc and thus appear to be 
less consistent with the observations by S04. The total 
particle numbers for coUisionlcss and coUisional parti- 
cles for the SMC are 100000 and 60000, respectively. 
The total particle number of a MC system is 260000 
(or larger, depending on whether the gas disk is mod- 
eled), when the LMC is also represented by a fully 
consistent disk model rather than a point-mass parti- 
cle. The total mass of the dark matter halo (Mdm) in a 
model with MsMC =4.5x10^Mq is about 2.4x IO^Mq: 
Mdm depends on /g and Msmc- 
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Figure 4: The top panel shows the time evolu- 
tion of the distances between the LMC-SMC (thick 
sohd), the LMC-Galaxy (thin solid), and the SMC- 
Galaxy (dotted) for the fiducial model Gl. The 
middle panel shows the time evolution of the to- 
tal masses within the central 5 kpc of the SMC 
for the dark matter (solid), old stars (dotted), gas 
(short-dash), and new stars (long-dash). The bot- 
tom panel shows the time evolution of the total 
masses the dark matter (solid), old stars (dotted), 
gas (short-dash), and new stars (long-dash) that 
are initially in the SMC and later in the central 
7.5 kpc of the LMC at a give time step. 



2.2.2 The disk model 

The radial (R) and vertical (Z) density profile of the 
initially thin disk of the SMC in this disk model arc 
assumed to be proportional to exp(— i?/i?o) with scale 
length Ro = 375 pc and to sech^(Z/Zo) with scale 
length Zo = 0.2Ro respectively. The stellar disk is em- 
bedded in a massive dark matter halo with the NFW 
radial density profile and the halo is truncated at the 
tidal radius of the SMC (n — 7.5 kpc). The mass ratio 
of the dark matter (Mdm) to the stellar disk (Msmc.s) 
is set to be 2.33 for all disk models. In addition to 
the rotational velocity made by the gravitational field 
of disk and halo component, the initial radial and az- 
imuthal velocity dispersion are given to the disk com- 
ponent according to the epicyclic theory with Toomre's 
parameter (Binney & Tremaine 1987) Q = 1.5. The 
vertical velocity dispersion at a given radius is set to 
be 0.5 times as large as the radial velocity dispersion 
at that point. 

The gaseous component in the disk model has the 
same properties as that in the spheroidal model. For 
consistency with the dwarf spheroidal models, we in- 
vestigate the disk models with 1 < /g < 2 and 2 < 
rg < 4. Thus, the initial structures and kinematics of 
old stars are quite different in the dwarf spheroidal and 
disk models for the SMC whereas those of the gas are 
the same in the two models. These initial differences in 
stellar components can result in significant differences 
in dynamical properties and star formation histories 
of the SMC between the two different models. The 
particle numbers for the three components in the disk 
models are the same as those in the dwarf spheroidal 
model. 

2.2.3 Star formation and chemical enrich- 
ment 

Star formation is modeled by converting the coUisional 
gas particles into collisionless new stellar particles ac- 
cording to the algorithm of star formation described 

below. We adopt the Schmidt law (Schmidt 1959) with 
exponent 7 = 1.5 (1.0 < 7 < 2.0, Kcnnicutt 1998) as 
the controlling parameter of the rate of star forma- 
tion. The amount of gas consumed by star formation 
for each gas particle in each time step is given as: 



P'g = C'sfPg'' 



(10) 



where pg and Csf is the gas density around each gas 
particle and a constant in the Schmidt law. The value 
of Caf is determined such that the star formation rate 
in the isolated model for the Galaxy with the disk mass 
of 6 X 1O^°M0 and the gas mass fraction of 0.1 can be 
~ IM0 yr"^ (e.g., Bekki 1998). We convert a gas 
particle into a field star only if the local surface gas 
density exceeds the observed threshold gas density of 
~ 3 Mq pc^^ for Magellanic dwarf irregular galaxies 
(Hunter ct al. 1998). We estimate the surface gas den- 
sity of a SPH particle from the volume gas density by 
assuming that the scale-height of the gas disk in the 
SMC being 100 pc during simulations. These stars 
formed from gas arc called "new stars" (or "young 
stars") whereas stars initially within a disk are called 
"old stars" throughout this paper. 
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Chemical enrichment through star formation and 
supernova feedback during the SMC evolution is as- 
sumed to proceed both locally and instantaneously in 
the present study. We assign the metallicity of origi- 
nal gas particle to the new stellar particle and increase 
the metals of the each neighbor gas particle with the 
total number of neighbor gas particles equal to A'gas, 
according to the following equation about the chemical 
enrichment: 

AMz = {ZiJinietms + (1.0-iimet)(1.0-Zi)msJ/met}/Afgas 

(11) 

here the AAfz represents the increase of metal for each 
gas particle. Zi, iimet, rns, and j/met in the above equa- 
tion represent the metallicity of the new stellar parti- 
cle (or that of original gas particle), and the fraction 
of gas returned to interstellar medium, the mass of the 
new star, and the chemical yield, respectively. The 
values of -Rmot, Umet, and the initial metallicity are 
set to be 0.3 and 0.004, and 0.0032, respectively. For 
these values, the final mean metallicity of of the SMC 
is [Fe/If] = —0.56: this value is significantly smaller 
for a smaller ymet. 

Dynamical evolution of interstellar medium (ISM) 
can be significantly influenced by supernova explosion 
and the effects of supernova explosion on ISM depend 
on the initial mass function (IMF) adopted in the present 
study. For the stellar population formed in a GMC, we 
assume an IMF described as i/'(mi) = Ami~'', where 
mi is the initial mass of each individual star and the 
slope of s = 2.35 corresponds to the Salpeter IMF. The 
normalization factor A is a function of md, which is a 
total mass of a stellar system, mi (lower mass cut-off), 
and mu (upper one): 



niu is a parameter and set to be lOOM© for most mod- 
els whereas mi is regarded as a key free parameter 
in the present study. We shows the results of the 
canonical models with mi = O.IMq and mu = lOOM© 
and the top-heavy IMF ones with mi = I.OMq and 
niu = IOOMq . The mass fraction of stars that can ex- 
plode as type-II supernova (i.e., rni > 8Mq) is 14% for 
the canonical model and 35% for the top heavy IMF 
ones. About 10 % of the total energy of one supernova 
(~ 10^^ erg) is assumed to be converted into kinematic 
energy of gas around the supernova (Thornton et al. 
1998). 

Some fraction of the kinematic supernova energy 
(EsN,i) of the ith stellar particle is assumed to be re- 
turned to gas particles (ISM) and then to change the 
kinetic energy of gas particles around the iih stellar 
particle. Each jih gas particle with the mass of mg j 
around an ith stellar particle can receive a velocity 
perturbation (Avsnj) directed radially away from the 
ith stellar particle. Avsnj satisfies the relation: 

fsNEsN,i = 2 X] "^9.j(^^SN,i)^, (13) 

i=i 

whore fsN represents the fraction of supernova energy 
returned to ISM, NndA is total number of gas parti- 
cles that are located within the smoothing {hi) radius 



around an ith stellar particle. Although we investigate 
the models with /sn = (no feedback), fsN ~ 0.1, and 
/sN = 1.0 (maximum feedback), we show the results of 
the models with the more realistic /sn = 0.1 (Thorn- 
ton et al. 1998). 

2.2.4 Orbital configurations 

The initial spin of a SMC (gas) disk in a model is spec- 
ified by two angles, 9 and </>, where 9 is the angle be- 
tween the Z-axis and the vector of the angular momen- 
tum of a disk and (f) is the azimuthal angle measured 
from X-axis to the projection of the angular momen- 
tum vector of a disk onto the X — Y plane. Although 
we investigate many models with different 9 and 4>, 
we mainly describe the results of the models with 9 = 
45° and (f) = 230° , This is mainly because these models 
can not only explain fundamental properties of the MS 
(GN and YN) but also reproduce dynamical properties 
of the present SMC. 

2.2.5 The GRAPE-SPH simulations 

We use the GRAPE (GRAvity pyPE) system (e.g., 
Sugimoto et al. 1991) for the smoothed particle hy- 
drodynamics (SPH) in order to evolve self-gravitating 
stellar and gaseous systems. Since we have already de- 
veloped our original "GRAPE-SPH" code capable of 
evolving self-gravitating systems composed both of col- 
lisionless and coUisional masses (Bekki &: Chiba 2006), 
we use the code to simulate the evolution of the MCs 
for the last ~ 3 Gyr. In this code, we calculate (i) 
mutual gravitational interactions between all particles 
and (ii) the list of neighbor particles for each SPH 
gas particle by using the GRAPE. All other calcula- 
tions related to hydrodynamics, star formation, chem- 
ical evolution, and orbital evolution of the MCs are 
done by host computers for the GRAPE. 

In order to implement SPH on the GRAPE sys- 
tem, we adopt numerical methods and techniques es- 
sentially the same as those used for the TREESPH 
method (e.g., Hernquist & Katz 1989, HK89). We 
adopt the spherically symmetric spline kernel proposed 
by Monaghan & Lattanzio (1985) and determine the 
smoothing lengths of SPH particles based on the num- 
bers of their neighbor particles (HK89). We use the 
same forms of artificial viscosity as those used in HK89 
(given in equations 2.22 - 2.25 of HK89). We adopt 
an isothermal equation of state for the gas and deter- 
mine the proper range of the initial gas temperatures 
(Ti) by running isolated disk models. We find that 
for the models with Ti = 10* K, which is reasonable 
for luminous disk galaxies like the Galaxy, the initially 
thin gas disks quickly become very thick owing to the 
higher gaseous pressure. A smaller Ti is necessary for 
the SMC with its smaller mass and, thus, weaker self- 
gravity in the vertical direction. Although we investi- 
gate the models with Ti = 250 K and 10^ K, we mostly 
describe the results of the models with Ti = 250 K. 

We adopt the predictor-corrector algorithm that is 
accurate to 2nd order in time and space in order to 
integrate the equations describing the time evolution 
of self-gravitating systems in the present GRAPE-SPH 
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code. The size of the time step (Sti) for each SPH 
particle is determined by the Courant condition for 
the particle. 6ti for a SPH particle is derived from 
its gas density, sound velocity, velocity gradient, and 
smoothing length according to the formula adopted in 
HK89. The minimum and maximum possible sizes of 
time stops for SPH particles are set to be 1.67 x 10^ yr 
and 6.85 x 10^ yr, respectively. We choose the above 
very small minimum possible size so that the time step 
size of every gas particle at a given time step can not 
be smaller than the minimum possible size: otherwise 
the simulation needs to stop. The size of the time step 
for a coUisionless particle (dark matter, old and new 
stars) is fixed at 6.85 x 10^ yr. 

The gravitational softening lengths (e) are differ- 
ent between coUisionless and coUisional particles and 

different between models with different initial SMC 
sizes (-RsMc)- For the fiducial model (described later), 
e = 1.70 kpc for coUisionless particles and e = 199 pc 
for coUisional (gaseous) ones. Since our main interest 
is on the Ikpc-scale structures and kinematics of the 
SMC, we consider that the resolution of the present 
simulation is high enough to investigate them. 



2.2.6 Stellar population synthesis models 

We derive the two-dimensional (2D) distributions of 

the projected B-band luminosity densities (Sl,b) of 
the SMC based on the spatial distributions and the 
age- and metallicity-distributions of old and new stars. 
In order to do this, we use the predictions of stellar 
population synthesis models by Vazdekis et al. (1996), 
which are available on his website. Since the observed 
metallicity of the stellar component of the SMC is with 
[Fe/H] = -0.73 ± 0.03 (Luck et al. 1998), we use the 
B-band luminosities and stellar mass-to- light-ratios for 
the SSP models with [Fe/H] ~ —0.68 and variously 
different ages in Vazdekis et al. (1996). We assume 
that the mean ages of old and new stars in the present 
simulation axe 5 Gyr and 0.5 Gyr, respectively, and 
thereby allocate the -B-band luminosity to each stellar 
particle. We chose this set of stellar ages so that we 
can show how the SMC looks on the sky if the SMC has 
different distributions of old and new stars. Since our 
main focus in this paper is structural and kinematical 
properties of the SMC, we do not intend to discuss 
how the appearance of the SMC on the sky depends 
on ages of old and new stars. 

We mainly investigate the projected 2D El,b dis- 
tribution of the present SMC for the region with O.lh < 
a < 1.9hand— 77° <S < —69° (corresponding roughly 
to the physical range of 4 kpc x 4 kpc) so that observa- 
tions can be compared with the simulations. We divide 
the region into 50 x 50 cells and thereby estimate Sl,b 
for each cell. We then smooth out the 2D El,b field by 
using the Gaussian kernel with the size 0.5 times the 
cell size so that wo can clearly see the global 2D Sl,b 
distribution without urnicccssarily dramatic cell-to-cell 
variations of El,b due to small numbers of particles in 
some cells. 



2.3 Main points of analysis 

In the present GRAPE-SPH simulations, the SMC is 
modeled as a fully self-gravitating N-body system whereas 
the LMC is represented by a test particle with the 
above-mentioned Plummer potential for most models. 
The mass distribution of the Galaxy is assumed to 
be fixed so that dynamical influences of the MCs on 
the Galaxy are not explicitly included in the models. 
We therefore do not intend to investigate dynamical 
and hydrodynamical interactions between the SMC, 
the LMC, and the Galaxy: these interactions will be 
investigated in detail by Bekki & Chiba (2008a), in 
which the influences of the LMC-SMC interaction on 
the LMC's evolution, and those of the LA-Galaxy in- 
teraction on the evolution of the outer part of the 
Galactic HI disk are discussed. 

We consider that the best model of the SMC in the 
present study needs to reproduce reasonably well the 
projected locations of the MS and the LA in the Galac- 
tic coordinate {l,b). Accordingly, we try to search for 
the models that can reproduce the observed locations 
of the MS and the LA among models with different 
model parameters. We first run "coUisionless mod- 
els" (labeled as "C") in which star formation and gas 
physics arc switched off. If we find that a coUisiorfloss 
model reproduces the locations of the MS and the LA 
well, we then run a "gas dynamical model" (labeled as 
"G") with model parameters exactly the same as the 
coUisionless model. Given that the present chemody- 
namical simulations are numerically costly, this two- 
fold investigation is quite effective in finding the best 
model of the SMC among G models. The LMC is rep- 
resented by a fully self-consistent disk model (paper 
I) in some coUisionless models so that we can con- 
firm that the formation of the MS does not depend 
on whether the LMC is modeled as a test particle or a 
self-consistent disk composed of many particles. 

Fig. 3 shows the distribution of "gas" particles 
with respect to the center of the Galaxy for the model 
CI, in which orbital evolution and masses of the MCs 
are consistent with previous models by GN and C06. 
This model can reproduce the observed location of the 
MS and the LA reasonably well and shows bifurcated 
structures in the MS, which appears to be consistent 
with observations. This model is therefore one of the 
successful models for the MS in the present study. The 
physical properties of the MS can be well reproduced 
by the present models with different Msmc , as long as 
we adopt the GN velocity type. However, the models 
with the K06 velocity type can not reproduce well the 
observed locations of the MS and the LAs, irrespective 
of Msmc- The physical properties of the MS and the 
LAs in models with different velocity types are briefly 
discussed in the Appendix A. 

Although we run 35 models with different model 
parameters, we describe the results for 14 represen- 
tative models that we regard as important in terms 
of the comparison with observations. The values of 
model parameters for 12 models are given in Table 2: 
The model numbers (column 1), AfsMC (2), velocity 
types (3), initial morphologies of the SMC (4), i?sMC 
(5), rs (6), /g (7), rg (8), Ti (9), and comments (10). 
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We mainly describe the results of the "fiducial model" 
in which the simulated properties of the SMC are more 
consistent with observations. Models with K06 veloc- 
ity types fail to explain the observed locations of the 
MS and the LA in the galactic coordinate system. We 
however show and discuss the results of these in order 
to discuss the latest proper motion measurements of 
the MCs (K06, Piatek et al. 2008) and orbital evolu- 
tion models of the MCs based on them (Besla et al. 
2007). 

We also investigate the "isolated model" with no 
LMC-SMC-Galaxy interaction in order to confirm that 
the derived characteristic properties of the SMC in 
the G models are due to the LMC-SMC-Galaxy in- 
teraction. We investigate three isolated models (II, 
12, and 13) with different Ti so that we can confirm 
which model(s) can be dynamically stable. We find 
that the gas disks in models with higher Ti (e.g., 13 
with Ti = 10* K) become unrealistically thick owing 
to high gaseous pressure not reasonable for the SMC 
with a small mass. We also find that the models with 
Ti = 250 K can keep initially thin gas disks and con- 
tinue to be stable. We therefore show the results of 
the G models with T — 250 K in the present study. 
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Figure 5: The final distributions of old stars (ma- 
genta) and new ones (cyan) on the sky (i.e., in the 
a — 6 coordinate system) for the fiducial model 
(Gl). 



3 Results 

3.1 Fiducial model 

3.1.1 Dynamical evolution 

Fig. 4 shows how the total masses of dark matter, old 
stars, new ones, and gas within the central 5 kpc of 
the SMC evolve with time during strong tidal inter- 
action between the SMC, the LMC, and the Galaxy 
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Figure 6: The final 2D distributions of projected 
mass densities (E in units of M© pc~^) for old stars 
(Es) and new ones (Ens) m logarithmic scales for 
the fiducial model. Note that the fiattencd struc- 
ture can be clearly seen in this simulated distribu- 
tions. 



in the fiducial model GI. The dark mater, old stel- 
lar, and gaseous components lose 40%, 41%, 75% of 
their initial masses during the tidal interaction. The 
larger fraction of gas mass lost during the interaction 
is due both to more efficient tidal stripping of gas and 
to gas consumption by star formation. The central 
spheroidal component (i.e., old stars) in this model can 
not be so strongly influenced by the tidal fields owing 
to its initial compactness so that it does not lose its 
mass during the strong LMC-SMC interaction around 
T = —1.5 Gyr. It however loses a significant fraction 
of old stars at the last LMC-SMC interaction about 
0.2 Gyr ago so that the MB composed of old stars can 
be formed. The final total mass within the central 5 
kpc in this model is 2.4 x lO^M©, which is consistent 
with the observed total mass of the SMC (e.g., S04). 

Fig. 4 also shows that not only stars but also gas 
initially in the SMC can pass through the central 7.5 
kpc of the LMC (corresponding to the disk size), in 
particular, at T = —1.4 and —0.2 Gyr. This result con- 
firms the "Magellanic squall" (Bekki & Chiba 2007), 
which means efficient mass-transfer from the SMC to 
the LMC via the LMC-SMC tidal interaction. Fig. 
4 furthermore shows that (i) new stars formed in the 
SMC can pass through the LMC and (ii) the total mass 
of the SMC's new stars that are located within the 
LMC at T = is about 4 x IO^Mq. These results 
imply that metal-poor ([Fe/H] < —0.6), younger stars 
(or star clusters) with ages less than 3 Gyr currently 
observed in the halo region of the LMC might well 
originate from the SMC. We note that Muhoz et al 
(2006) recently discovered metal-poor stars with pos- 
sible metallicities of [Fe/H] = —0.67 in the outer halo 
of the LMC, though they did not estimate ages of these 
stars. Owing to the total mass of the LMC being sig- 
nificantly larger than that of the SMC, mass transfer 
from the LMC to the SMC does not occur. 

The strong tidal interaction can also compress strongly 
the gas to form high-density asymmetric structures 
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Figure 7: The final projected mass density profiles 
(S) in units of Mq deg~^ along a for old stars 
(solid), new ones (dotted), and gas (dashed) in 
the fiducial model. The profiles are derived for 
for particles with a > asMC) where asMC is a 
of the SMC, so that the mass profiles in the MB 
regions can be clearly seen. For comparison, the 
MB region where HI and molecular observations 
(e.g., Mizuno et al. 2006) have been done is shown 
by two dotted lines. 



(e.g., spirals), wfiero star formation is moderately en- 
hanced (up to ~ O.4M0 yr~^). Owing to the absence 
of the stellar bar in the dwarf spheroid, a nuclear star- 
burst is not induced by the interaction: this result is 
in a striking contrast with those by Noguciii (1988) in 
which stellar bars formed in the disks of interacting 
galaxies can trigger nuclear starbursts. About 30% of 
the initial gas is tidally stripped to form the MS, the 
LAs, and the MB, depending on the present locations 
of the gas. The stellar spheroid can not be strongly 
distorted by the tidal interaction so that the final mor- 
phology of the SMC shows an old stellar spheroid with 
a gaseous disk more compact than the initial one and 
a compact stellar one composed only of young stars 
with ages less than 2.7 Gyr. 

3.1.2 Stellar structures and kinematics 

As shown in Fig. 5, the final projected spatial distribu- 
tion on the sky (i.e., on the a-5 coordinate system) are 
quite different between old and new stars. New stars 
originate from the gas disk disturbed and compressed 
by the tidal fields of the LMC and the Galaxy so that 
the final distribution of new stars is quite elongated: 
this structure is reminiscent of the observed apparently 
"barred structure" in the SMC. The old stars, on the 
other hand, arc not so strongly influenced by the tidal 
fields as the gas owing to their initial compact distri- 
bution so that they can not be so tidally elongated at 
r = 0. They appear to show an almost spherical shape 
on the sky, though the spherical morphology is not so 
clearly seen in this plot because of the stripped stars 
being widely distributed in the outer halo region of the 
SMC. These result clearly explain why there are signif- 
icant differences between the morphology of old stel- 
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Figure 8: The same as Fig. 7 but along S. The 
profiles are derived for for particles with S > Ssmc, 
where 5smc is 5 of the SMC, so that the observed 
mass profiles of the stellar halo can be compared 
with the simulation. The physical ranges of a and 
5 in Figs. 7 and 8 are the same so that the profiles 
along a and 5 can be compared with each other. 



lar populations (e.g., rod giant stars; Harris & Zarisky 
2006) and the optical morphology depending strongly 
on young stellar populations. A GC-like compact stel- 
lar system that is formed during the tidal interaction 
and thus composed of new stars can be seen in the up- 
per left corner of the plot (i.e., (q, 5) « (1.5h, —70.5°)). 

Fig. 6 shows that the 2D distributions of the pro- 
jected surface mass densities (E) for old stars (Es) and 
for new ones (Ens) that are produced based on the 
projected distributions of stars for the region shown 
in Fig. 5. The old stellar component clearly shows 
a spherical distribution on the sky, though Eg of the 
outer halo is slightly enhanced in the northern direc- 
tion. This north-south asymmetry in the outer stellar 
halo is highly likely to be due to the past SMC-LMC- 
Galaxy tidal interaction. The new stellar component 
is significantly elongated and the major axis of the 
distribution in the central region (_R < 1 kpc) is not 
coincident with that of the outer part {R > 2 kpc). 
The lower part appears to be more strongly disturbed 
than the upper one owing to the last LMC-SMC tidal 
interaction about 0.2 Gyr ago. 

Owing to the efficient tidal stripping of old stars 
in the last LMC-SMC interaction about 0.2 Gyr ago, 
the SMC loses a significant fraction (~ 40%) of old 
stars. As a result of this, the stripped stars can be 
widely distributed between the LMC and the SMC 
(i.e., the MB region). Fig. 7 shows that the pro- 
jected surface mass density (E) of old stars in the MB 
region has a relatively flat distribution with E being 
significantly higher than lO^'Modeg^'^. The presence 
of this extended halo-like component is due largely to 
the stripped. The surface density of new stars in the 
MB region is significantly lower than that of old ones 
in this model, though the E-ratio of old to new stars 
in the MB depends on the initial distributions of old 
stars in the SMC, as discussed later in this paper. As 
shown in Fig. 7, there are no significant differences in 
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Figure 9: The final projected 2D distribution of 
the i?-band luminosity density (Sl,b) of the SMC 
for the fiducial model. The distribution is pro- 
duced based on the mass distributions shown in 
Fig. 6 for adopted stellar population models. The 
details of the methods to produce the distribution 
are given in the main text. 
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the E profiles between gas and new stars in the MB. 

Fig. 8 shows that the E profiles of old and new 
stars and gas along 5 appear to be significantly dif- 
ferent from those along a, which suggests that the 
global mass distribution in the outer part of the SMC 
is highly asymmetric owing to the past interaction 
with the LMC and the Galaxy. E of gas and new 
stars along a are systematically higher that those along 
5 for the outer halo region (|a — asMcl > 0.7h or 
\S — 5sMc| > 3deg), which means that a larger frac- 
tion of the stripped stars and gas are located in the 
MB region. There are no remarkable differences in E 
profiles along a and 5 for old stars in the SMC. These 
results on the possible differences in E profiles along a 
and 5 between stellar and gaseous components will be 
discussed later in the context of the formation of the 
MB and the outer stellar halo in the SMC. 

Fig. 9 shows that the 2D B-band surface bright- 
ness (El,b) distribution of the SMC has an irregular 
shape rather than a regular spherical one, which is ex- 
pected from the dominant old stellar populations in 
this model. This result clearly suggests that the ob- 
served irregular appearance in the SMC is due largely 
to the distribution of young stars formed during the 
tidal interaction with the LMC and the Galaxy. How- 
ever, the simulated B-band morphology is not so sim- 
ilar to the observed optical morphology, which has 
a "wing" and a non-elliptical appearance in its cen- 
tral region. These apparent differences between the 
simulated and observed morphologies imply that the 
present models are not so successful in reproducing the 



Figure 10: The final profiles of line-of-sight veloc- 
ities {Vr — Vr,sMC: solid) and dispersions (cr, dot- 
ted) for old stars (top), new ones (middle), and 
gas (bottom) in the fiducial model, where Vr and 
K-.SMC are the line-of-sight velocity of particles and 
of the center of the SMC, respectively. For com- 
parison, the observed range of the velocity disper- 
sions of older stars is shown by dashed lines. 



2D distributions of star-forming regions in the SMC. 

Fig. 10 shows that the line-of-sight velocity dis- 
persion (a) profile for old stars has a strong radial 
gradient with the maximum a of 37 km s~^ and the 
mean a of 31 km s~^. The line-of-sight velocity pro- 
file {Vr — Vr,SMc) for old stars appears to have a small 
amount of rotation {V < 20 km s~^). These simulated 
kinematical properties of the old stars are broadly con- 
sistent with the observed ones, a is systematically 
lower in new stars than in old ones, which reflects 
that the new stars originate from gas that can dis- 
sipate away the random kinematical energy produced 
during the SMC-LMC-Galaxy tidal interaction. New 
stars have a significant amount of rotation {V ~ 40 km 
s~^) and thus a high V/a (> 1). These results suggest 
that young stellar populations with ages less than ~ 3 
Gyr in the SMC can have kinematical properties signif- 
icantly different from older ones such as PNe, carbon 
stars, and red giants. 

As shown in Fig. 10, the line-of-sight velocity and 
velocity dispersion profiles of the gas are clearly differ- 
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Figure 11: The AMR of new stars formed for 
the last - 2.7 Gyr evolution of the SMC in the 
fiducial model. The la dispersions in metallicities 
are shown by error bars for each of the age bins. 




Figure 13: The same as Fig. 6 but for the disk 
model G3 in which old stars are modeled as a ro- 
tating disk. 



-0.2 

X 



-0.4 



i - 



X -0.6 - 

td 

a 

-0.8 h-,-;-,-T-,-r-,-,-|-,-;-,-;-rrr,-;-r-;-,-,--^i 

1 2 3 4 5 

Radius (kpc) 

Figure 12: The radial dependence of the mean 
metallicities of new stars in the SMC for the fidu- 
cial model. The la dispersions in metallicities are 
shown by error bars for each of the radial bins. 



ent from those of the old stars in the sense that the 
gas show lower velocity dispersions (7 — 42 km s~^), 
a higher maximum rotational velocity (~ 49 km s^^), 
and thus a higher V/a of ~ 1.2. The derived kine- 
matical difference between old stars and gas and the 
high maximum rotational velocity are broadly consis- 
tent with the latest kinematical observations by S04. 
Fig. 10 also shows that the velocity profiles for stellar 
and gaseous components are significantly asymmetric 
in the sense that the absolute magnitudes of the ve- 
locities appear to be systematically larger for a > Ih. 
This asymmetry might well be due partly to the last 
LMC-SMC interaction about 0.2 Gyr ago. 

3.1.3 Stellar populations 

In this model with a relatively high j/met = 0.004, 
the mean metallicity of the new stars evolves from 
[Fe/H] = -0.75 at T = -2.2 Gyr to [Fe/H] = -0.56 
at T = Gyr (note that the initial gaseous metallic- 
ity is [Fe/H] = —0.8). It should be stressed here that 
the final metallicity of new stars can be significantly 



smaller than the above value for smaller i/met: the ob- 
served mean metallicity of young stars with ages less 
than ~ 3 Gyr can give constraints on the chemody- 
namical models of the SMC. As shown in Fig. 11, there 
is an age-metallicity-relation (AMR) of new stars (i.e., 
more metal-rich for younger stars), though the disper- 
sions in stellar metallicities appears to be quite large, 
especially for younger stars. The shape of the AMR 
and the metallicity dispersions also strongly depend on 
j/mct, as discussed later in this paper. 

Fig. 12 shows that the mean metallicities of new 
stars are higher in the inner regions than in the outer 
ones for the central 3 kpc of the simulated SMC. This 
negative metallicity gradient of new stars reflects the 
fact that chemically enriched gas due to star forma- 
tion is transfered to the central region and converted 
into new stars there during the SMC-LMC-Galaxy in- 
teraction. The metallicity dispersion is largest in the 
innermost radial bin, which implies that new stars cur- 
rently in the central region originate from gas that 
are initially located in different regions with different 
metallicities. Although the negative metallicity gradi- 
ent of new stars is a robust prediction of the present 
models, the slope of the gradient and the metallicity 
dispersions at radial bins depend strongly on i/mct. 

3.2 Parameter dependences 

Parameter dependences of stellar and gaseous proper- 
ties are summarized as follows. 

(1) Final spatial distributions of stars do not de- 
pend strongly on model parameters in the dwarf spheroidal 
models: old and new stars have spherical and flat- 
tened distributions on the sky, respectively. The de- 
rived structural differences between old and new stars 
are broadly consistent with observations, which implies 
that the SMC was born as a dwarf spheroidal. 

(2) The flnal structural and kinematical properties 
are quite different between the dwarf spheroidal and 
disk models. Fig. 13 shows the final spatial distri- 
butions of old and new stars projected onto the sky 
for the disk model G3. The spatial distribution of old 
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Figure 14: The same as Fig. 9 but for the model 
G3. 



stars do not appear to be so spherical, which suggests 
that dynamical heating by the tidal fields of the LMC 
and the Galaxy is not strong enough to morphologi- 
cally transform the initially disky shape into the more 
spherical one. New stars in the disk model show a 
bar in the projected distribution on the sky, which 
is consistent with the observed optical morphology of 
the SMC. Although only old stars appear to form the 
outer extended stellar halo, the spatial distributions 
between old and new stars are quite similar with each 
other. Fig. 14 shows that the final B— band morphol- 
ogy in the disk model can be identified simply as an 
inclined disk galaxy. 

(3) The final velocity dispersions (a) of old stars 
are roughly similar {a ~ 30 km s~^) between different 
dwarf spheroidal models. The maximum line-of-sight- 
velocities (V) of old stars are less than 20 km s~^ for 
most of the dwarf spheroidal models. It should be 
here stressed that the apparent non-negligible amount 
of rotation {V < 20 km s~^) in old stars are due to the 
streaming motion of the stars in the outer region of the 
SMC. The final velocity dispersions and maximum ro- 
tational velocities of new stars are systematically lower 
and higher, respectively, than those of the old ones, ir- 
respective of model parameters. 

(4) There are no significant kinematical differences 
between old stars, new ones, and gas in the disk mod- 
els, which are inconsistent with observations. Fig. 
15 shows that these three components show higher V 
(~ 40 km s~^) in the disk model G3 even after the 
strong tidal heating of the SMC by the LMC and the 
Galaxy. The velocity dispersions of new stars are only 
slightly smaller than those of old stars in these disk 
models. These results suggest that the disk models 
are less consistent with observations than the dwarf 
spheroidal ones in terms of kinematical properties of 
the SMC. 
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Figure 15: The same as Fig. 10 but for the model 
G3. 



(5) The surface mass densities (E) in the MB region 
depend strongly on initial distributions of stars of the 
SMC in the present models. Fig. 16 shows an example 
of models in which E of old stars can be significantly 
lower than those of new ones and gas. The model G2 
shows that E of old stars is lower than lO^M© deg~^ 
in some regions of the MB. This is mainly because a 
smaller number of old stars are tidally stripped in this 
model owing to the more compact initial distribution 
of the stars. These results imply that the observed E 
of stars in the MB region can give strong constraints 
to the models for the initial distribution of old stars in 
the SMC. 

(6) The SMC can lose a significant fraction (up to 
~ 40% depending on model parameters) of its initial 
mass in the present models. As a result of this, the 
final masses of the SMC can be significantly smaller 
than the observationally suggested one 2.4x lO^Mo 
within the central 3 kpc of the SMC; S04) for some 
models (e.g., G6, CI, and C3) with Msmc = 3.0 x 
10^ Mq: the initial sizes of the SMC (-Rsmc) need 
to be quite small so that too efficient tidal stripping 
can be prevented in these low-mass models. Fig. 17 
shows that in the compact, low-mass model, the SMC 
loses only 26% of its initial mass so that the final 
mass within the central 5 kpc (= 2.1 x lO^M©) can 
be similar to the observed mass of the SMC (S04). 
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Figure 16: The same as Fig. 7 but for the model 
G2 in which the initial spheroidal component is 
more compact than that in Gl. 



The present models with Msmc = 3.0 x lO^M© and 
^SMC = 7.5 kpc shows significantly lower final masses 
(< 2.0 X lO^M©), which suggests that the previous MS 
models in C06 may have serious problems in reproduc- 
ing the present SMC mass. Fig. 17 confirms that the 
Magellanic squall (i.e., gas transfer to the LMC from 
the SMC during the LMC-SMC interaction) can be 
seen in this low-mass SMC model. 

(7) The final mean stellar metallicity, the AMR of 
new stars, and the radial stellar metallicity gradient in 
the SMC depend on j/mot. For example, as shown in 
Fig. 18, the AMR in the model G4 with ymct = 0.002 
is significantly flatter than that in Gl. Both the stel- 
lar metallicity for young stars with ages less than 0.2 
Gyr and the dispersions in stellar metallicities in this 
model are significantly smaller than those in Gl. It is, 
however, currently unclear which of the two models, 
Gl and G4, is more consistent with observations on 
abundances of young stars with ages less than 3 Gyrs. 

4 Discussion 

4.1 Origin of stellar and gaseous kine- 
matics in the SMC 

Previous kinematical studies based on radial velocities 
of older stellar populations such as PNe and red giants 
suggested that the stellar component of the SMC has 
the velocity dispersion of 20 — 30 km s"'^ and appears 
to have no global rotation (e.g., Dopita et al. 1995; 
Suntzeff et al. 1986; Hatzidimitriou et al. 1993). The 
latest survey of 2046 red giant stars have suggested 
that the stellar component of the SMC is primarily 
supported by its velocity dispersion (Harris & Zaristky 
2006). Observational studies of HI gas, on the other 
hand, showed that the gaseous component has global 
rotation with the peak rotational velocity of about 60 
km s~^ (S04). No theoretical studies, however, have 
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Figure 17: The same as Fig. 4 but for the model 
G5 in which both the initial total mass and the 
size of the SMC are significantly smaller than that 
in Gl. 



yet provided a physical explanation for this intrigu- 
ing kinematical difi^erence between stellar and gaseous 
components in the SMC. 

The present study has first shown that kinemati- 
cal differences between old stars and gas can be clearly 
seen only in the dwarf spheroidal models: both stel- 
lar and gaseous components of the disk models show 
significant rotation even after strong tidal heating by 
the LMC and the Galaxy and have no kinematical dif- 
ferences between them. These results imply that the 
SMC previously did not have a rotating disk composed 
both of older stars and gas but had an extended gas 
disk kinematically decoupled from its stellar compo- 
nent. The rotating gaseous disk is one of essential 
ingredients in the previous successful models for the 
MS (GN, YN, and C06). Therefore the formation pro- 
cesses of the MS described in these studies would not 
be changed significantly if stellar spheroids rather than 
stellar disks were adopted in these studies. The previ- 
ous coUisionless models (BC07b) and the present ones 
in which the SMC is assumed to be a dwarf spheroidal 
can reproduce physical properties of the MS reason- 
ably well (see Appendix A for more discussions on the 
MS formation). 

The present study suggests that younger stars formed 
from the rotating gas disk of the SMC for the last ~ 2.7 
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Figure 18: The same as Fig. 11 but for the model 
G4 with ymet = 0.002. 



Gyr can sliow a smaller (projected) central velocity 
dispersion {a ^ 20 km s~^) and a larger (projected) 
rotational velocity {V ~ 40 km s~^). Although 3D 
structures of Cophoid variables (possible candidates of 
younger stars) along the line-of-sight were previously 
investigated for the SMC (e.g., Caldwell & Coulson 
1986), kinematical properties of young stellar popular 
tions are largely unknown. It is thus doubtless worth- 
while for future observations to confirm whether or not 
younger stellar populations with ages less than 3 Gyrs 
in the SMC show rotational kinematics unlike its older 
stellar populations. 

4.2 Formation of a dwcirf spheroidal 
with an extended HI disk 

The present study has first shown that the observed 
dynamical properties of the SMC can be successfully 
reproduced by the new dwarf spheroidal models. The 
previous disk models which have been often adopted 
in explaining the MS (e.g., YN03) would not be rea- 
sonable. We therefore suggest that the SMC before in- 
teracting strongly with the LMC and the Galaxy (i.e., 
more than 3 Gyr ago) was a gas-rich dwarf like the 
NGC 6822 with both a spheroidal component com- 
posed of old and intermediate-age stars and an ex- 
tended HI disk (e.g., de Blok & Walter 2003). If the 
SMC was really a gas-rich dwarf spheroidal, then the 
following question arises: how were the main body of 
the central spheroid and the outer HI gas disk formed 
? 

Mayer et al. (2001) showed that low-luminosity 
disk galaxies in the Local Group can be transformed 
into dwarf spheroidal and elliptical galaxies, if they are 
strongly influenced by tidal fields of the giant galaxies 
like the Galaxy. Since this transformation scenario 
predicts quite effective tidal stripping of gas and stars 
in the outer parts of disks, it would not be possible 
that the transformed dwarf spheroids can still contain 
the outer HI disk. Therefore the formation of the SMC 
with an extended HI disk would not result from past 
tidal interaction of the SMC with the Galaxy. 

Recent HI observations of the SMC have suggested 
that the total mass of the SMC derived from the HI ro- 



tation curve is about 2.4 x 10"Mq and thus consistent 
with a dynamical model without a dark matter halo 
(e.g., S04). These observations may well imply that 
the SMC is a "tidal dwarf" which can be formed from 
tidal tails during strong tidal interaction between two 
galaxies (e.g.. Due et al. 2000), because tidal dwarfs 
can not contain dark matter components that are ini- 
tially in halos of their host galaxies. Although there 
are two luminous galaxies (i.e., M31 and the Galaxy) 
in the Local Group which can be host galaxies for tidal 
dwarfs, there appears to be no observational evidence 
for their past strong interaction. Thus the tidal dwarf 
scenario currently appears to be less promising as a 
formation model of the SMC. 

Thus we suggest an "accretion scenario" in which 
the SMC is likely to have developed its spheroidal 
component at the epoch of galaxy formation and then 
gradually formed its extended HI disk via gaseous ac- 
cretion from its outer halo. The potential problem of 
this accretion scenario is that the strong tidal field of 
the LMC and the Galaxy could not have allowed the 
halo gas of the SMC to be accreted onto the SMC's 
disk, if the SMC has been bound closely to the LMC 
and the Galaxy since its formation. The present tidal 
radius of the present SMC is only ~ 5.0 — 7.5 kpc (GN, 
YN, and C06), which implies that only a small amount 
of gas could have been accreted onto the SMC with- 
out being stripped by the LMC and the Galaxy. The 
present distances of the SMC with respect to the LMC 
and the Galaxy may have been significantly larger in 
the past than in the present one. If this is the case, the 
SMC may have developed its massive HI disk through 
gradual accretion from the halo to finally become a 
dwarf spheroidal with an extended HI disk. Other sce- 
narios (e.g., major merging) for the formation of the 
SMC are discussed in detail by Bekki & Chiba (2008b). 

In the above discussions, we assume that the SMC 
has a dynamically hot spheroid with no/little rotation. 
It is however possible that the SMC is actually a rota- 
tionally flattened spheroid viewed almost from face-on 
(i.e., along the rotational axis): we can not detect ro- 
tation of the spheroid in the line-of-sight-velocity field, 
even if the old stars of the SMC have rotational kine- 
matics. If this is the case, the SMC is almost like a 
polar-ring galaxy in which the rotational axes of stars 
and gas are perpendicular with each other. Although 
this possibility can not be ruled out by current ob- 
servations, the probability of such a configuration of 
the SMC with respect to the Galaxy would be very 
low. Thus we suggest that the observed apparent dif- 
ferences in kinematics between old stars and HI gas in 
the SMC reflect real differences in internal kinematics 
between them: such differences may well suggest that 
the SMC would have experienced accretion/merging 
event of very gas-rich objects (e.g., gas-rich dwarfs) in 
its dynamical history possibly in a small group where 
galaxy merging is more likely. 

4.3 On the possible absence of older 
stars in the MB 

Observational studies revealed a number of intriguing 
physical properties of stars and gas in the inter-Cloud 
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regions and the MB (e.g., Irwin ot aL 1990; Mullcr 
et al. 2004; Mizuno et al. 2006; Harris 2007; Kato 
at al. 2007; Nishiyama et al. 2007). Muller et al. 
(2004) reported that the MB is composed of two kine- 
matically and morphologically distinct structures and 
shows a remarkable velocity offset around 6 = —73.5° 
in the S — Vrad piano. The discovery of molecular clouds 
(Mizuno ct al. 2006) and the possible detection of 
prc-rriain-sequence stars (Nishiyama ct al. 2007) in 
the MB strongly suggest that star formation in the 
MB is currently ongoing and possibly triggered by the 
last LMC-SMC interaction. Harris (2007) did not de- 
tect an older stellar population in the MB and thus 
suggested that the MB originates purely from gaseous 
components of the MCs. Since the observed structural 
and kinematical properties of the MB and the possi- 
ble physical mechanisms for star formation in the MB 
have been already discussed by Muller & Bekki (2007) 
and Bekki & Chiba (2007), respectively, we here focus 
on the origin of the possible absence of older stellar 
populations in the MB. 

The present model with the GN velocity type and 
Msmc,s ~ IO^Mq predicts that the projected surface 
mass density in the MB is roughly 6 x 10* — 3.0 x 
lO^'Modcg^"^. The presence of old stars in the MB is 
due largely to the tidal stripping of stars during the 
last LMC-SMC close encounter with its periceriter dis- 
tance of 8 kpc about 0.2 Gyr ago. The presence of 
old stars in the MB thus appears to be inconsistent 
with observations by Harris (2007), though the obser- 
vational upper limit of the mass density of an older 
stellar population is not clearly shown in Harris (2007). 
The following two are the possible interpretations for 
this apparent inconsistency between the simulations 
and the observations. The first is that the simulated 
last LMC-SMC interaction needs to be weaker, so that 
fewer old stars end up in the MB. The second is that 
the initial size of the stellar spheroid (or disk) in the 
simulation needs to be smaller to prevent its efficient 
tidal stripping by the LMC. 

Noel & Gallart (2007) have recently found the pres- 
ence of intermediate-age and old stars in the southern 
regions at 6.5 kpc from the SMC center. This result, 
which is in a striking contrast to the results by Har- 
ris (2007), suggests that the SMC has am outer stel- 
lar halo component and is possibly more massive than 
other observations have ever suggested. Furthermore, 
this result possibly implies that there could be older 
stars in the MB, if the stellar halo of the SMC has 
a spherical distribution. The upper limit for the pro- 
jected mass density of older stars in the MB is crucial 
to determine whether the last LMC-SMC interaction 
is the main mechanism of the MB formation. We thus 
suggest that more extensive observational studies need 
to be done to provide strong constraints for the forma- 
tion models of the MB. 

If the projected surface density of older stars in the 
MB is truly low, as the result by Harris (2007) sug- 
gests, the result by Noel & Gallart (2007) possibly im- 
plies that the stellar halo of the SMC is quite asymmet- 
ric and elongated to the north-south direction. The 
strong tidal fields of the LMC and the Galaxy can 
transform the outer shape of the SMC's stellar halo 



so that the recent interaction history, which is deter- 
mined by the orbits of the MCs, can be imprinted on 
the shape of the stellar halo. We thus suggest that 
the direction of the major axis of the possibly elon- 
gated stellar halo might well provide some constraints 
on the orbit of the SMC around the LMC: the direc- 
tion would be perpendicular to the orbit of the SMC 
and point to the center of the LMC (thus can constrain 
the orbital plane of the LMC-SMC system). 

4.4 Dark matter properties of the 
SMC 

We have first pointed out that (i) the dark matter halo 
of the SMC is highly likely to have a very large core 
(> 5 kpc) and (ii) the observed rotation curve of HI 
and the total mass of the SMC within the central few 
kpc can be explained better by the SB profile of the 
dark matter distribution. The SB profile with a large 
core might well explain why the SMC appears to have 
no dark matter in the central few kpc (S04), though 
a fully self-consistent dynamical model including the 
dark matter halo has not been construed yet for the 
SMC (Bekki & Chiba 2008a). If the SMC has a such 
large, low-density core of its dark matter halo, it is 
highly Ukely that the halo was initially quite extended 
and therefore lost a substantial amount of its original 
mass owing to the past tidal interaction with the LMC 
and the Galaxy since its formation. 

As shown in the present simulations, significant 
fractions of original masses in the SMC (up to ~ 40%) 
can be lost during the last 2.7 Gyr LMC-SMC-Galaxy 
interaction. We have therefore suggested that in or- 
der to explain the observed total mass of the SMC 
within the central 3 kpc, the original mass of the SMC 
should be significantly more massive than Msmc = 
3.0 X W^Mq, which is adopted in previous MS mod- 
els (GN and COG). We have demonstrated that the 
models with Msmc = 4.5 x lO^M© can better explain 
the present SMC mass even after the significant mass 
loss of the SMC due to the interaction. It should be 
noted here that such models with larger Msmc and 
thus larger HI gas mass can naturally solve the prob- 
lems of the previous MS models in explaining the ob- 
served total mass of the MS and the LAs (e.g., YN03). 
We thus suggest that the SMC initially has a larger 
total mass (possibly as much as 4.5 x lO^M©) and a 
large core in the dark matter distribution. 

5 Conclusions 

We have investigated structural, kinematical, and chem- 
ical properties of stars and gas in the SMC interact- 
ing with the LMC and the Galaxy based on chemody- 
namical sinmlations. We have adopted a new dwarf 
spheroidal model with an extended HI gas disk for 
the SMC and investigated the chemical and dynam- 
ical evolution and the star formation history of the 
SMC for the last ~ 2.7 Gyr. The main results are 
summarized as follows. 

(1) The final spatial distribution of the old stars 
projected onto the sky is more spherical even after the 
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strong LMC-SMC-Galaxy interaction whereas that of 
the new stars is significantly flattened and appears to 
be a bar. This difference in spatial distributions be- 
tween the old and new stars is due to the fact that 
new stars are formed in the thin HI disk of the SMC 
during the tidal interaction. These result suggest that 
the observed "bar" is mostly composed of now stars re- 
cently formed in the SMC. Some fraction of old stars 
initially in the SMC are stripped to form the outer 
stellar halo of the SMC. The distribution of the stel- 
lar halo projected onto the sky appears to be slightly 
asymmetric in the north-south direction, which possi- 
bly reflects the strong tidal interaction with the LMC 
about 0.2 Gyr ago. 

(2) Old stars have a line-of-sight velocity dispersion 
(cr) of ~ 30 km s~^ and weak rotation with the max- 
imum rotational velocity (V^) of less than 20 km s"'^ 
and show asymmetry in the profiles. New stars have 
a smaller a than old ones and a significant amount 
of rotation (V/a > 1). These kinematical differences 
between old and new stars do not depend on model 
parameters. It should be stressed that the derived V 
of old stars are highly likely to be overestimated by 
the stream motions of the stars along the tidal tails of 
the SMC. 

(3) The HI gas shows velocity dispersions of cr ~ 
10 — 40 km s^^ depending on the distances from the 
SMC's center, a high maximum rotational velocity {V ~ 
50 km s~^), and a spatial distribution similar to that 
of new stars. The derived higher velocity dispersion of 
gas is due to the tidal disturbance of gas by the LMC 
and the Galaxy. 

(4) The mean metallicity of stars can increase from 

[Fe/H] 0.80 to [Fe/H] 0.56 for y^^t = 0.004. 

The stellar populations show a weak AMR in the sense 
that younger stars are more likely to be more metal- 
rich. The final steUar [Fe/H] and the AMR depend 
strongly on j/met (e.g., the final metallicity is [Fe/H] ^ 
—0.66 for t/mot ~ 0.002). More metal-rich stars are 
more likely to be located in the inner regions of the 
SMC, irrespective of model parameters. 

(5) Although the MB can be formed in most mod- 
els, the simulated MB inevitably contains old stars. 
This is mainly because some fraction of old stars ini- 
tially in the spheroidal component of the SMC are 
tidally stripped to be distributed in the MB region. 
The mean surface density of the stars in the MB is 
different in different models and ranges from ~ 6 x 
lO^Mgdeg"^ to ~ 3 X lO'^Mgdeg"^. If these numbers 
are inconsistent with observations, these results mean 
that the orbital evolution of the MCs for the last 0.2 
Gyr in our models needs to be significantly modified so 
as to reduce significantly the stellar mass in the MB. 

(6) The disk model shows flattened spatial dis- 
tributions and higher maximum rotational velocities 
{V > 40 km s^^) in old stars, which appear to bo in- 
consistent with observations. Thus the dwarf spheroidal 
model can better explain observations than the disk 
one for a given orbit of the SMC in terms of structures 
and kinematics of old stars in the SMC. Chemical en- 
richment processes are not so different between the 
dwarf and the disk models. 

(7) The simulated MS, which is composed mostly 



of HI gas, shows bifurcation in its structure and kine- 
matics. The models based on the GN velocity type 
can much better explain distributions and kinematics 
of HI gas in the MB than those based on the K06 one 
that are derived from the latest HST proper motion 
measurements. Since we did not investigate the MS 
formation for all of the possible unbound orbits of the 
MCs, it is unclear whether all of the unbound orbits 
for the MCs consistent with K06 are highly unlikely to 
explain the MS formation. These results imply that 
the physical interpretations of the HST proper motion 
measurements need to be reconsidered in the context 
of the MS formation. We have suggested that future 
abundance studios of the MS through QSO absorption 
lines can provide valuable information on formation 
processes and early chemical evolution of the outer HI 
disk of the SMC. 

(8) The SMC is suggested to initially have a large 
core (> 5 kpc) in its dark matter distribution and a 
larger total mass (Msmc > 3.0 x IO^Mq) so that the 
observed kinematical properties of HI in the present 
SMC can be better explained. The more massive SMC 
with a low-density dark matter halo is suggested to 
explain better observations of the inner profile of the 
HI rotation curve and of the total gas mass in the MS 
and the LAs. 

The present study adopted the "classical" orbits of 
the MCs (GN) rather than the new one (K06) in order 
to compare results of the present new SMC model with 
those of the previous ones (e.g., Yoshizawa & Noguchi 
2003; C06) that also adopted the classical orbits. We 
thus have not solved the problems related to the incon- 
sistency between the predicted velocities of the MCs 
and the observed ones (K06). In our forthcoming pa- 
pers, we discuss how to solve the problems based on 
more sophisticated numerical simulations. 
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A Formation of the MS 

The previous tidal models of the MS reproduced suc- 
cessfully a number of observed physical properties of 
the MS, such as the presence of the LA, the absence 
of stars in the MS, and bifurcated structures and kine- 
matics of the MS (e.g., MF, GN, YN, and C06). Al- 
though the ram pressure models, in which the MS orig- 
inates from the gas disk of the LMC, provided some 
explanations for the origin of the observed kincmatical 
properties of the MS, they failed to explain the pres- 
ence of the LA and the bifurcated structure of the MS 
(e.g., Mastropietro et al. 2005). Given that the sim- 
ple version of the tidal models has some difficulties in 
explaining in the location of the LA on the sky (Bokki 
et al. 2008), it would be fair to claim that it remains 
unclear whether the tidal interaction between the MCs 
and the Galaxy is the dominant formation mechanism 
for the MS. Also it is unclear whether the present new 
dwarf spheroidal models can explain physical proper- 
ties of the MS as well as the previous ones. 

Recent observational studies based on proper mo- 
tion measurements of the MCs by HST ACS have shed 
new light on the above unresolved problem of the MS 
formation (K06, Besla et al. 2007). Based on a sys- 
tematical study of orbital evolution of the MCs, Besla 
et al. (2007) suggested that the MCs arrived at the 
Galaxy about 0.2 Gyr ago for the first time in their 
histories, if the observed proper motions of the MCs 
are correct. These studies thus imply that the MS is 
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Figure 19: The final distributions of gas in 

the galactic coordinate system (Lb) for the dwarf 
spheroidal model CI. This model is designed to 
compare with the early results obtained by the col- 
lisionless MS models of GN and C06 in which the 
SMC is modeled as a disk galaxy. 



unlikely to be formed as a result of the past tidal in- 
teraction of the SMC with the LMC and the Galaxy 
about 1.5 Gyr ago. It is, however, unclear whether 
theoretical models with the orbits of the MCs being 
consistent with observations by K06 can reproduce the 
observed properties of the MS. 

Thus, we need to investigate both (i) whether the 
observed MS can bo well reproduced by the present 
new SMC models with the GN velocity type and (ii) 
whether the present models with the K06 one can ex- 
plain the observed properties of the MS as well. We 
show the results of the two collisionless models (CI and 
C2), because other models show similar properties to 
those of these two models. The detailed results of the 
MS formation in more self-consistent N-body models 
for the MCs and the Galaxy will be described in Bekki 
& Chiba (2008a) 

Fig. 19 shows that the locations of the MS and the 
LAs in the model CI are broadly consistent with the 
observed ones (e.g., Putman et al. 1998), though the 
shapes of the LAs appears to be too broadened. The 
bifurcation of the MS can be clearly seen, which con- 
firm the early results by C06. Fig. 20 clearly shows the 
velocity gradient of the MS along the Magellanic lon- 
gitudes, which is consistent with observations. Fig. 20 
also shows the bifurcation in kinematics (i.e., Vlsr) for 
the Magellanic longitude larger than 50 degrees. These 
simulated bifurcation in the structural and kinematical 
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Figure 20: Vlsr (local-standard-of-rest velocities) 
of gaseous particles as a function of Magellanic lon- 
gitude for the model CI. 



properties of the MS will soon be compared with the 
corresponding observations (e.g., Briins et al. 2005; 
Nidever et al. 2007) so that the validity of the models 
will be addressed. 

Fig. 21 shows that the simulated distributions of 
gas on the sky are inconsistent with observation in 
the model C2 with the K06 velocity type, though two 
gaseous streams can be formed from the extended gas 
disk of the SMC as a result of SMC-Galaxy interac- 
tion. It is also confirmed in the present study that the 
observed locations of the MS and the LA arc not re- 
produced by the present models with the K06 velocity 
types irrespective of model parameters (e.g., Msmc). 
There can be mainly two interpretations of the above 
results. The first is that the although the 3D space 
motions of the MCs derived from the observed proper 
motions of the MCs arc very close to the true ones, the 
adopted assumption that the MS originates from the 
SMC is wrong. Besla et al. (2007) proposed a scenario 
in which the MS originates from gaseous components 
initially within the LMC possibly through gaseous out- 
flow triggered by stellar feedback effects. Owing to the 
lack of numerical studies based on the above scenario, 
it remains unclear whether the MS can originate from 
the LMC, if the 3D space motions of the MCs proposed 
by K06 are correct. 

The second is that the proposed 3D space motions 
of the MCs by K06 arc not close to the true ones so that 
the inconsistency in the 3D space motions between the 
tidal models (e.g., GN) and the observationally sug- 
gested ones (K06) can not be a serious problem for 
the tidal models. Bekki (2008) first showed that the 
mean proper motion of the SMC derived from the three 
local fields can be significantly different from that of 
the center of the mass for the SMC so that the 3D 
proper motion of the SMC in K06 can deviate signif- 
icantly from the true one: interpretations of the ob- 
served proper motion measurements may or may not 
be correct. These results imply that the second inter- 
pretation is reasonable, though it remains still unclear 
how close the 3D space motions adopted in the tidal 
models are to the true ones. 
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Figure 21: The same as Fig. 19 but for the model 
C2. 



Whether the first or the second interpretation is 
correct, a number of physical properties of the MS have 
not been clearly explained. In particular, the "kink" 
in the LA around 6 = 0° (C06) and the "Magellanic 
filaments" (Nidever et al. 2007) have not been repro- 
duced by any numerical models of the MS formation 
in a self-consistent manner. As suggested by Bekki et 
al. (2008), the origin of the kink around b = 0° can re- 
sult from the hydrodynamical interaction between the 
outer HI disk of the Galaxy and the LA. The apparent 
physical connections between the Magellanic filaments 
and the LMC imply that there could be some interac- 
tions between the MS and the LMC. We plan to discuss 
these observations in Bekki & Chiba (2008a) based 
on more sophisticated models in which both a LMC- 
MS (or LMC-LA) hydrodynamical interaction and a 
Galaxy-LA one can be self-consistently investigated. 

If the MS originates from the SMC, the gaseous 
abundances should be very similar to those of the outer 
gas disk of the SMC about 1.5 Gyr ago. The present 
chemodynamical models imply that the gaseous metal- 
licity in [Fe/H] ranges from —0.8 to —0.6, which is 
significantly smaller than that of the LMC ([Fe/H] ~ 
—0.3, van den Bergh 2000). The details of abundance 
patterns (e.g., [C/Fe], [N/Fe], and [Mg/Fe]) of the MS 
should be similar to those of the SMC about 1.5 Gyr 
ago. We thus suggest that future extensive observa- 
tional studies of abundance patterns of gas in the MS 
from QSO absorption lines will provide vital clues to 
the question as to whether the MS originate from the 
SMC or from the LMC. 



